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ABSTRACT 

We use numerical simulations of structure formation in a Cold Dark Mat- 
ter (CDM) cosmology to compare the angular momentum distributions of dark 
matter and non-radiative gas in a large sample of halos. We show that the two 
components have identical spin parameter distributions and that their angular 
momentum distributions within individual halos are very similar, all in excellent 
agreement with standard assumptions. Despite these similarities, however, we 
find that the angular momentum vectors of the gas and dark matter are poorly 
aligned, with a median misalignment angle of ~ 30 degrees, which might have 
important implications for spin correlation statistics used in weak lensing stud- 
ies. We present distributions for the component of the angular momentum that is 
aligned with the total angular momentum of each halo, and find that for between 
5 and 50 percent of the mass this component is negative. This disagrees with 
the generally adopted 'Universal' angular momentum distribution, for which the 
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mass fraction with negative specific angular momentum is zero. We discuss the 
imphcations of our results for the formation of disk galaxies. Since galactic disks 
generally do not contain counter-rotating stars or gas, disk formation cannot oc- 
cur under detailed conservation of specific angular momentum. We suggest that 
the material with negative specific angular momentum combines with positive 
angular momentum material to build a bulge component, and show that in such 
a scenario the remaining material can form a disk with a density distribution 
that is very close to exponential. 

Subject headings: cosmology: dark matter — galaxies: formation — galaxies: 
structure — galaxies: halos. 



1. Introduction 

Understanding the structure and formation of disk galaxies is intimately linked to under- 
standing the origin of its angular momentum. In hierarchical structure formation scenarios 
the luminous parts of galaxies form from gas that is cooling and condensing within dark 

matter (DM) halos and these merge to build larger and larger objects (White & Rccs 1978). 
Within this picture, the current paradigm for disk formation contains three important ingre- 
dients: (i) the angular momentum originates from cosmological torques (Hoyle 1953), (ii) the 
gas and dark matter within virializcd systems have initial angular momentum distributions 
that are identical (Fall & Efstathiou 1980), and (iii) the gas conserves its specific angular 
momentum when cooling (Mestel 1963). 

It is well established that cosmological torques impart angular momentum to dark mat- 
ter halos. Many studies have investigated the resulting angular momenta of dark matter 
halos using either tidal torque theory (e.g., Peebles 1969; Doroshkevich 1970; White 1984; 
Catelan & Theuns 1996) or N-body simulations (e.g., Barnes & Efstathiou 1987; Warren et 
al. 1992; Bullock et al. 2001). The specific angular momentum of halos is typically parame- 
terized by the dimensionless spin parameter 

with G the gravitational constant and J, E, and M the total angular momentum, energy and 
mass of the object, respectively (Peebles 1969). Under the assumptions (n) and (iii) disks 
have scale lengths Rd oc Ait!vir- Given the distributions of halo spin parameters and of halo 
virial radii i?vir, the implied distribution of disk scale lengths is in excellent agreement with 
observations (e.g., Dalcanton, Spergel & Summers 1997; Mo, Mao & White 1998; de Jong 
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& Lacey 2000). This success has prompted many detailed studies of disk galaxy formation, 
always under the three assumptions listed above (van den Bosch 1998, 2000, 2001, 2002; 
Jimenez et al. 1998; Natarajan 1999; Heavens & Jimenez 1999; Firmani & Avila-Reese 2000; 
Avila- Reese & Firmani 2000; Buchalter, Jimenez & Kamionkowski 2001). 

The standard picture of disk formation that has emerged from these studies has been 
remarkably successful in explaining a wide variety of observational properties of disk galaxies. 
However, two significant problems, both related to the angular momentum, have come to 
light. First of all, detailed hydro-dynamical simulations of disk formation in a cold dark 
matter (CDM) Universe yield disks that arc an order of magnitude too small (Navarro 
& Benz 1991; Navarro & White 1994; Steinmetz & Navarro 1999; Navarro & Steinmetz 
2000). This problem, known as the angular momentum catastrophe, is a consequence of the 
hierarchical formation of galaxies which causes the baryons to lose a large fraction of their 
angular momentum to the dark matter. It is unclear whether this is a problem for the theory 
or merely for the particular way in which feedback processes arc implemented (or ignored) in 
the simulations (e.g., Navarro & White 1994; Weil, Eke & Efstathiou 1998; Sommer-Larsen, 
Gelato & Vedel 1999; Thacker & Couchman 2001). 

The second problem concerns the actual density structure of galaxy disks. Assuming 
detailed angular momentum conservation, this structure is a direct reflection of the distri- 
bution of specific angular momentum in the proto-galaxy. Bullock et al. (2001) determined 
the angular momentum distributions of individual dark matter halos, which according to 
assumption (ii) should also reflect the angular momentum distribution of the gas that forms 
the disk. However, these distributions seem to have far too much low angular momentum 
material to be consistent with the typical exponential density distributions of disk galaxies 
(Bullock et al. 2001; van den Bosch 2001; van den Bosch, Burkert & Swaters 2001). Knebe 
et al. (2002) and Bullock, Kravtsov & Colm (2002) investigated the angular momentum 
profiles of warm dark matter (WDM) halos. Except for a slight difference in the average 
spin parameter, CDM and WDM halos have very similar angular momentum distributions, 
implying that this particular problem cannot be solved by invoking WDM instead of CDM. 
A more promising solution has recently been suggested by Mailer & Dekel (2002). In their 
picture halo angular momentum originates from the orbital angular momentum of merging 
satellites, rather than from pure tidal torques (see also Mailer, Dekel & Somerville 2001). 
Combining this picture with a simple model for feedback and tidal stripping Mailer & Dekel 
are able to explain the angular momentum profiles of dwarf galaxies measured by van den 
Bosch et al. (2001). 

These angular momentum problems seem to imply that assumptions (ii) and (iii) listed 
above, and which have been essential ingredients of our standard paradigm for the formation 
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of disk galaxies, may be incorrect. Remarkably, despite many papers that have investigated 
the angular momentum distributions within dark matter halos, to the best of our knowledge, 
no study so far has addressed the spin distribution of the gas in these halos in order to test 
whether the ansatz that the gas and dark matter have similar initial distributions of specific 
angular momentum is correct. In fact, it seems plausible that they may not. Gas and 
dark matter suffer different relaxation mechanisms during halo collapse: Whereas the dark 
matter undergoes coUisionless virialization through violent relaxation, the gas attempts to 
achieve a hydrostatic equilibrium through shocks. In order to test whether the initial angular 
momentum distributions of gas and dark matter are indeed similar, we perform numerical 
simulations of structure formation in a ACDM cosmology including both dark matter and 
a non-radiative gas. After presenting our simulations and outlining our analysis methods 
(Section 2), we present various statistics of the angular momentum distributions (Section 3). 
In Section 4 we discuss the implications of our results for the theory of disk galaxy formation, 
and we summarize our results in Section 5. 



2. Method 

2.1. Numerical Simulations 

We use the public version of the smoothed particle hydrodynamics (SPH) code GAD- 
GET developed by V. Springel (1999) and described in Springel, Yoshida & White (2000). 
The chosen flat cosmology is vacuum energy dominated with JIa = 0.7 and matter content 
Q,jn — 0.3 with a baryon density Q;, = 0.021/i~^ and a Hubble constant of 67kms~^ Mpc~^ 
(i.e., h — 0.67). We choose a relatively small box of 10/i~^ comoving Mpc and evolve it 
from rcdshift 59 down to redshift three. We use identical power-spectra for baryons and 
dark matter given by the fitting formula of Efstathiou, Bond & White (1992) and use 128^ 
particles each for the gas and the dark matter component. The power spectrum is normal- 
ized to have a typical mean overdensity of (jg = 0.9 at rcdshift zero in a top hat sphere of 
8/i~^ Mpc. Gas and dark matter particles have masses of mgas = 6.26 x IQ^'h^^ M© and 
m-oM — 3.34 X 10^/i~^ M©, respectively. The gravitational softening lengths for both the gas 
and the dark matter were chosen to be Ah~^ kpc (comoving). The gas has an initial entropy 
corresponding to a temperature of 2 x 10^ K for a mean molecular weight of one proton 
mass. An artificial temperature fioor at lO'^ K (for the same molecular weight) is introduced 
to mimic the heating effects from re-ionization. Both this temperature fioor and the initial 
temperature are too low to play a significant role in the simulated dynamics because they 
are significantly lower than the virial temperatures of the halos we analyze. This is evident 
in our analysis as the reported results do not depend on the mass of the halos studied. All 
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results presented below correspond to the final output at 2; = 3. 

The robustness of our results with respect to resolution, simulation box size, and the 
redshift of analysis is addressed in Section 3.3. 

2.2. The Halo Identification Algorithm 

We start by using the publicly available group finder HOP, written by Eisenstein & Hut 
(1998), to identify overdensity peaks in the dark matter distribution using a threshold density 
of 20 times the mean dark matter density. In each group we identify the densest particle, 
and we determine the virial radius it!vir of the spherical volume, centered on this densest 
particle, inside of which the average dark matter density is Avii.(2;) times the critical density 
PcTiti^)- For our adopted cosmology and at the redshift of the output analyzed Ayir = 174.7 
(Bryan & Norman 1998). Next we compute the center-of-mass of all dark matter particles 
inside this spherical volume which we associate with the new center of the halo. We compute 
the new i?vir around this new halo center, as well as the associated new center-of-mass. This 
procedure is iterated until the distance between the center-of-mass and the adopted halo 
center is less than one percent of the virial radius i?vir- Typically this requires of the order 
of 2 to 5 iteration steps. 

All dark matter and gas particles inside the resulting spherical volume with radius -Rvir 
are considered halo members. Although our iterative method assures that the center-of- 
mass of the dark matter component is similar to the adopted halo center, the same is not 
necessarily true for the gas component. We therefore remove all halos from our sample for 
which the distance between the halo center and the center-of-mass of the gas particles is 
larger than 10 percent of i?vir- We checked that none of our results correlate with the offset 
between the centers of mass of the gas and dark matter. In addition, we repeated the entire 
analysis adopting the position of the most bound particle as halo center and found very 
simular results. 

In order to allow sufficiently accurate measurements of the angular momentum vectors 
of the gas and dark matter we only accept halos that have more than 100 gas particles and 
more than 100 dark matter particles. In addition, if any two halos overlap, i.e., the distance 
between the halo centers is less than the sum of their virial radii, we remove the least massive 
halo from our sample. This leaves a total of 378 halos of which 259 are isolated in that they 
have no subgroups detected, and they do not overlap with any other halo. Of these 378 halos 
23 have more than 5000 particles (in total) , and are used for a more detailed analysis of the 
distribution of specific angular momentum (see Section 3.2). 
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2.3. Theoretical Background 



Our main goal is to compare the angular momentum distributions of the gas and dark 
matter components. However, one cannot do so by simply computing angular momentum 
distributions from the velocities of the gas and dark matter particles. Both the dark matter 
and the gas are fluids, for which the velocity of each microscopic particle can be written 
as V = u + w. Here u is the mean streaming motion at the location x of the microscopic 
particle, and w is the particle's random motion. The macroscopic, coUisionless dark matter 
particles in i¥-body simulations can be considered a Monte Carlo realization of the phase- 
space distribution function of all microscopic particles. The velocities of these particles in the 
simulation correspond to v and are computed by solving Newton's equations of motion. In 
the case of the gas, however, the Smoothed Particle Hydrodynamics (SPH) approach is used, 
which yields the streaming motions u of the gas particles by solving the Euler equations (with 
some artificial viscosity to take account of shocks). Information about the random motion 
of the gas particles is provided by the internal energy of each particle. 

Thus the velocities of the dark matter particles and gas particles in the simulation 
correspond to different motions. This must be taken into account before comparing their 
angular momentum distributions. This is most easily achieved by computing the velocities 
V of the gas particles from their streaming motions and their (isotropic) velocity dispersion 
tensor. For each gas particle we randomly draw 100 random velocities w which we add to 
the particle's streaming motion u. For each of the three Cartesian components of w we draw 
a velocity from a Gaussian with a standard deviation given by 



Here U and T are the internal energy per unit mass and temperature of the gas particle, 
respectively, k is Boltzmann's constant, and n is the mean molecular weight of the gas. After 
this 'thermal broadening' we find that the rms velocities (|vp) of the gas and dark matter 
particles are, for the vast majority of our halos, very similar, indicating that the halos are, 
to a good approximation, in hydrostatic equilibrium. 

We thus use the velocities v to compare the angular momentum distributions of the gas 
and dark matter. However, the distribution that is of interest for galaxy formation is the 
angular momentum distribution based on the streaming motions u(x). After all, when the 
gas cools V — > u, i.e., the typical random motions become negligible, at least in sufficiently 
massive halos. Under the assumption that the gas conserves its specific angular momentum, 
the resulting disk will thus have an angular momentum distribution related to its streaming 
motions u. Ideally one would therefore like to compare the angular momentum distributions 
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computed from the streaming motions u(x). However, this requires estimating the streaming 
motions of the dark matter. In principle this could be achieved by smoothing the dark 
matter velocities v, but it is unclear what smoothing scale to use and resolution issues are 
likely to play an important role. We therefore do not attempt to compute the dark matter 
streaming motions, but instead use the velocities v of the dark matter and gas particles 
when comparing their angular momentum distributions. However, the streaming motions 
of the gas are available, and we present a detailed comparison of the angular momentum 
distributions of the gas based on v and u in Section 3.2. 



2.4. Analysis 

For each halo we compute the total angular momenta of the gas and the dark matter 
particles, using, 

-^^gas.DM 

Jgas,DM = "T'gas.DM ^ X (Vj — Vcom)- (3) 

i=l 

Here Vcom is the mean (bulk) velocity of the entire halo (dark matter plus gas) and is the 
radial vector with respect to the halo center. We follow Bullock et al. (2001) and define the 
following modified spin parameters 

^/ _ |jgas,DM| / ,s 

V2i?virKi/ ^ ^ 

Here jgas and Jdm are the specific angular momenta of the gas and dark matter, respectively, 
and Kir = a/ ^(Mgas + Mdm) / -Rvir IS the circular velocity at the virial radius i?vir- This 
definition of the spin parameter has the advantage that it does not depend on the energy 
content of the halo which can be difficult to estimate reliably. In what follows we always 
adopt the spin parameter definition of equation (4), and we drop the prime. 

In addition to the spin parameters of the dark matter and the gas we also compute the 
angle 

(5) 

between the angular momentum vectors of the two components. After rotating the coordinate 
frame such that the 2;-axis is aligned with the angular momentum vector of either the dark 
matter or the gas, we also compute the fractions of dark matter and gas particles with 
negative angular momentum about their respective z-directions. These fractions are denoted 
by /dm and /gas, respectively, and give an indication of the amount of disordered motion 
relative to the ordered motion. 



e = cos"^ 



Jgas ■ Jdm 



'gas I 



"dm 
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Note that since at each position the expectation value of w is, by definition, the null 
vector, the values of Agas, Adm and 9 are independent of whether one considers the streaming 
motions u or the microscopic velocities v. However, the angular momentum distributions 
may be quite different. In order to distinguish these two cases we label the appropriate 
quantities with the superscript v whenever the microscopic velocities v are used. Thus a 
particle's specific angular momentum is given by either j = r x u or j'' = r x v. 



Figure 1 plots the histograms of Adm and Agas for all 378 halos in our sample. Both 
distributions are well fit by a log-normal distribution 



For the dark matter we find Adm = 0.040 and ct^dm — 0-56, while for the gas Agas — 0.039 
and — 0.57. Clearly, the distributions of Agas and Adm are indistinguishable. If we only 
consider the 259 isolated halos we obtain very similar distributions both for the dark matter 
and for the gas. 

In the upper left panel of Figure 2 we plot Agas versus Adm- Although Agas and Adm 
are well correlated, with a Spearman rank coefficient = 0.67, which has a probability of 
virtually zero to occur under the null hypothesis of no correlation, the scatter is relatively 
large. The average of Agas — Adm is —0.001 with a standard deviation of 0.020, which confirms 
once again that on average the spin parameters of the gas and dark matter are very similar. 
However, on a one-to-one basis the two components can have quite different specific angular 
momenta; i.e., the average (median) of |Agas — Adm| is 0.015 (0.011). We find no significant 
trend between Agas — Adm and Myir (see upper left panel of Figure 4) indicating that the 
scatter is real, and not an artifact of discreteness noise. 

The upper right panel of Figure 2 plots /|as versus f^y^. With {fZu/ fgJ = 0.97 ± 0.10 
the two mass components have very similar mass fractions with negative specific angu- 
lar momentum. The middle panels of Figure 2 show that the values of /gag and /dm are 
strongly anti-correlated with their respective spin parameters. A similar anti-correlation is 
also present for the gas when /gas, is used instead of /gag (lower left panel). Again, none 
of these results change if we only focus on the 259 isolated halos. These anti-correlations 
indicate that it is the amount of ordered motion that is directly responsible for the amplitude 
of the total specific angular momentum. 

Finally, the lower right panel of Figure 2 plots /gas versus f^^. For all halos /gas < /gas, 
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Results 
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illustrating the trivial result that the streaming motions of the gas are more ordered than 
their microscopic motions. Note, however, that there is not simply a constant offset between 
/gas and Rather, the difference between /gas and f^^ is larger for halos with larger 

Agas, and the distribution of /gas is much broader than that of f^^. All these results have 
implications for the formation of disk galaxies, as we discuss in more detail in Section 4. 

3.1. Misalignment of Angular Momentum Vectors 

Figure 3 plots the distribution of the misalignment angle 9 (in degrees) between the 
total angular momentum vectors of the dark matter and the gas. The mean (~ 36°) and 
median (~ 27°) of the distribution are similar for the entire sample and for the subsample 
of isolated halos. 

Contrary to what one might expect, there is no correlation between the misalignment 
angle 6 and the absolute difference between Agas and Adm- This is evident from the upper 
right panel of Figure 4. The solid lines indicate the mean and the 68 percent interval of 
9. We do, however, find a small trend that the average misalignment angle increases with 
decreasing halo mass and spin parameter (see lower panels of Figure 4). We stress that there 
is no correlation between Myir and Agas so that these are two independent trends. Both are 
likely due to discreteness effects, in that the direction of the angular momentum vectors 
can be more accurately determined if either the total angular momentum or the number of 
halo particles is larger. However, since both trends are weak and the average misalignment 
angle is still significantly different from zero even for the most massive halos and for the ones 
with the largest spin parameters, we conclude that there is a true misalignment between the 
angular momenta of the gas and the dark matter (e.g., see column (9) of Table 1). 

This misalignment may have important implications for galaxy spin correlation statis- 
tics. In the past, both tidal torque theory (Crittenden et al. 2001; Heavens, Refegier & 
Heymans 2000; Natarajan et al. 2001) and pure N-body simulations (Croft & Metzler 2000) 
have been used to calculate the possibility of intrinsic alignments of galaxy shapes, which 
may be responsible for part of the measured weak lensing signal (e.g. van Waerbecke et 
al. 2000; McKay et al. 2001), and may be used to reconstruct the shear field (Lee & Pen 
2000, 2001). In these studies the assumption is made that the angular momentum vectors 
of the galaxies are ahgned with those of the dark matter halos. The relatively large mis- 
alignment angles found here, however, may diminish any alignment of galaxy shapes that 
one might expect based on tidal torque theory. We emphasize, however, that simulations 
of higher resolution are required to reduce the possible impact of discreteness effects on the 
distribution of misalignment angles. 
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3.2. The Distribution of Angular Momentum 



Although ~ 100 particles is enough to obtain a reliable estimate of the spin parameter, 
many more particles are required to determine angular momentum distributions (hereafter 
AMDs) within individual haloes. We therefore limit ourselves here to halos with more than 
5000 particles in total, of which there are 23 in our sample. 

We define the normalized specific angular momentum / = j/(-RvirK'ir) and compute 
separate distributions p{l) for dark matter particles and for gas particles. Here and in what 
follows j refers to the specific angular momentum along the 2;-axis. Note that i?virKir is the 
maximum specific angular momentum that a halo particle can have (i.e., I e [—1, 1]). We 
normalize p{l) to unity so that 



In Figure 5 we plot p{l^) both for the gas (dashed fines) and for the dark matter (solid 
fines) of 20 halos that are selected randomly from the sample of 23. Some global parameters of 
these 20 halos are listed in Table 1. As can be seen, the AMDs of the two mass components 
are remarkably similar. The fact that the p(/^) of the gas appear less noisy than for the 
dark matter is a consequence of the fact that because of the 'thermal broadening' of the 
gas velocity field we efi^ectively have on average 100 times more gas particles per halo than 
dark matter particles (see Section 2.3). In a few cases small differences can be discerned, but 
there is no indication for any systematic difference between the AMDs of the two components. 
Therefore, the standard assumption that the gas and dark matter have the same distributions 
of specific angular momentum is to first order correct, and the shocks that occur during the 
virialization process do not significantly modify the angular momentum distribution of the 
gas. 

In Figure 6 we compare the gas AMD p{l'') (dashed lines) to the distribution p{l) as 
computed from the streaming motions of the gas (solid lines). Here a clear systematic 
difference is apparent; typically p{l) has a less extended wing to negative specific angular 
momentum than p{l'") (cf. lower right panel of Figure 2), and a more pronounced peak at 
low angular momentum. Note that with equation (7) both distributions integrate to the 
same value of Agas. 

In Figures 7 and 8 we plot X{R) and A9{R) as functions of R/Ryir, respectively. Results 
for both the dark matter (solid lines) and the gas (dashed lines) are shown. Here X{R) is 
the spin parameter inside a sphere of radius R, i.e.. 




(7) 



y/2GRM{R) 



(8) 
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with j'"{R) and M{R) the specific angular momentum and mass inside a sphere of radius 
R, respectively. A6{R) is defined as the angle between the angular momentum vector inside 
radius R and that inside i?vir (i-e., by definition A^(i?vir) = 0). In some cases the gas and 
dark matter have fairly similar X{R) and A9{R) profiles (e.g., halos 14, 19, and 58), while in 
other cases the gas and dark matter seem to have quite distinct profiles (e.g., halos 36, 156, 
582). Remarkably, halos for which the X{R) and/or A6{R) profiles are quite different for 
the gas and dark matter are not necessarily also the ones that have the least similar AMDs, 
nor do they correspond to halos for which any of the global parameters listed in Table 1 are 
in any way special. Note that A typically either decreases or stays roughly constant with 
radius, while the angular momentum directions of both the gas and the dark matter typically 
change quite dramatically from the center to the virial radius. We discuss the implications 
of this intrinsic misalignment for the formation of disk galaxies in Section 4. 



3.3. Robustness of results 

The results presented above are based on a numerical simulation with relatively low 
resolution (128^ particles), with a small box size (Lbox = 10h~^ Mpc comoving), and which 
has only been run to redshift z = 3. Too small resolution causes shocks to be poorly resolved, 
which could result in an underestimate of the possible angular momentum redistribution of 
the gas during its relaxation process. A small box size may cause an inaccurate representation 
of the cosmological tidal field responsible for the angular momentum of halos. Finally, it is 
not a priori clear that angular momentum distributions at 2; = 3 are representative for those 
at 2; = 0. 

In order to address the robustness of our results against these effects we analyze a 
high resolution numerical simulation of a cluster kindly provided by Naoki Yoshida. In 
this simulation exactly the same cosmological parameters are adopted as in the simulation 
discussed above (see Section 2.1). It is based on the technique of 'zooming in' on a region of 
interest (Tormen, Bouchet & White 1997). The parent simulation from which this region is 
chosen is the GIF-ACDM simulation of Kauffmann et al. (1999) which followed 256^ particles 
within a comoving box with Lbox = 141.3/i~^ Mpc. The Lagrangian region that collapses to 
form the second most massive halo in this GIF simulation (with virial mass ~ 8 x 10^^h~^MQ) 
has been resimulated with greatly increased mass and force resolution and including a non- 
radiative gas. This high resolution region is represented by ~ 2.2 x 10"'' dark matter and gas 
particles each, with masses of 1.2 X lO^^h-^ Mq and 2.1 x lO^/i^^ M©, respectively. Outside 
this high-resolution region ~ 3.1 x 10^ particles with larger masses are used to represent the 
cosmological tidal field. More details regarding this simulation can be found in Springel et 
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al. (2001) and Yoshida et al. (2000a,b) 

We analyze the angular momentum distributions of the dark matter and gas of this high 
resolution halo at 2; = 0. Using the standard friends-of-friends algorithm (Davis et al. 1985) 
we identify the most massive group and its most bound particle. Next the virial radius 
i?vir of the spherical volume, centered on this particle, is determined adopting Avir = 101 
(appropriate for this cosmology and rcdshift; Bryan & Norman 1998). The halo thus defined 
(hereafter referred to as halo 'SO') contains ~ 1.2 x 10^ particles (more than twice that 
of the biggest halo analyzed in the previous section), and has a virial mass and radius of 
Mvir = 8.7 X 10^'^h~^ Mq and 1.86h'^ Mpc, respectively. Analyzing the angular momentum 
distributions of the dark matter and gas particles as in Section 2.4 we obtain the spin 
parameters and mass fractions with negative specific angular momentum listed in Table 1. 
The resulting angular momentum distributions p{l'") and p{l), as well as X{R) and A9{R), 
are shown in Figure 9. 

Note that the angular momentum distributions of halo 'SO' are in no way systematically 
different from those of the sample of halos analyzed above. In particular, once again we 
find that the gas and dark matter have virtually identical angular momentum distributions 
p(P), and that a large mass fraction (~ 45 percent) of the gas has negative specific angular 
momentum. The A(i?) and A6{R) profiles are similar to that of halo 11 (cf. Figures 7 
and 8). Although one has to be careful to draw conclusions based on a single high resolution 
halo, these results suggest that the characteristics of the angular momentum distributions 
presented in Section 3.2 are not influenced by the relatively low resolution, by the small box 
size, or by the fact that we analyzed the halos at z = 3. 



3.4. Do Dark Matter Halos follow a Universal Angular Momentum Profile? 

Bullock et al. (2001; hereafter BOl) computed the angular momentum distributions of 
dark matter halos and argued that these follow a universal profile given by: 

p(.) = •^") ^ KJ/jms^) (g. 

Here M(< j) is the halo mass with specific angular momentum less than j, jmax is the 
maximum specific angular momentum of the entire AMD, and is a free fitting parameter. 
Note that P{j) corresponds to the cumulative distribution of the specific angular momentum. 
In order to compute these AMDs, BOl split up their halos in a number of spatial cells. The 
angular momentum profiles P{j) are then constructed by ranking the cells according to their 
total specific angular momentum j from which the best-fitting value of is determined. Note 
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that this cell averaging should in principle average out the random motions w of the dark 
matter particles, and equation (9) should thus be interpreted as the cumulative distribution 
of p{j) rather than p{j^). 

However, there are two problems. First of all, it is unhkely that the particular cell 
geometry adopted by BOl yields the true streaming motions of the dark matter (see the 
discussion in Section 2.3). This means that, depending on the actual cell geometry adopted, 
the AMD of equation (9) does not necessarily reflect the true p{j). Secondly, because j 
corresponds to a projected component the total j of any cell can be negative. BOl ignored 
this and removed all cells with negative j from their ranked list. Wc find that the p{l) of our 
halos are well fit by equation (9) for / > 0, in agreement with BOl. However, whereas BOl 
have implicitely assumed that there is no mass with negative specific angular momentum, 
we find large mass fractions with Z < 0, at least for the gas. Since measuring the true 
streaming motions of the dark matter is non-trivial we have not attempted to obtain the 
p{l) for the dark matter. However, since the p{l^) of the gas and dark matter arc virtually 
identical, it seems reasonable to assume that the same applies for the p(/)-distributions^. If 
confirmed, between ~ 5 and 50 percent of the dark matter mass has negative specific angular 
momentum, in clear disagreement with the 'Universal' profile of BOl. Furthermore, since the 
mass fraction with negative specific angular momentum is strongly anti-correlated with the 
total specific angular momentum, the 'Universal' angular momentum profile of equation (9) 
lacks an essential characteristic of the true AMDs. 



4. Implications of the Formation of Disk Galcixies 

In the standard picture of disk formation, (i) the specific angular momentum of disk 
galaxies arises from cosmological torques, (ii) is identical to that of their dark matter halos, 
and (iii) is conserved during the process of disk formation. Recently several authors have 
pointed out that this picture is inconsistent with observations, as it would result in disks that 
are too centrally concentrated (BOl; van den Bosch 2001; van den Bosch et al. 2001). This 
is a reflection of the fact that the specific angular momentum distributions of dark matter 
halos, as described by the 'Universal' profile of BOl, contain relatively too much low angular 
momentum material. 

As we argued in Section 1 this seems to suggest that one or more of the standard 
assumptions is incorrect. Our results clearly show that in the non-radiative case the angular 



•^Thc accuracy of this assumption needs to be checked using simulations of higher resolution than the one 
analyzed here. 
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momentum distributions of the gas and dark matter are virtually identical, thereby lending 
strong support to (ii). However, our results also indicate that assumption (iii) cannot be 
correct. As shown, the angular momentum distributions p{l) of the gas contain relatively 
large mass fractions with negative specific angular momentum. Disk galaxies, on the other 
hand, generally do not contain material with negative specific angular momentum. This 
suggests that detailed conservation of specific angular momentum has to be violated if one 
wants to build a disk out of the typical AMDs of the gas in proto-galaxies. 

Once the gas starts to cool, its density increases, enhancing the frequency of shocks 
through which the gas parcels exchange angular momentum. Gas parcels with negative jz 
can share their angular momentum with parcels with positive jz, thus increasing the mass 
fraction with low angular momentum. This, however, can only aggravate the low angular 
momentum problem. Apparently real disks do not contain a large fraction of low angular 
momentum material (van den Bosch et al. 2001). It is tempting, however, to speculate that 
this low angular momentum material might actually build a (near-zero angular momentum) 
bulge component (cf. van den Bosch 1998; Kepner 1999). The mass of such a bulge depends 
on the details of the angular momentum exchange during cooling. The mass with negative 
specific angular momentum may combine with very little mass but with high positive specific 
angular momentum (yielding a minimum bulge mass), or with a relatively large amount of 
mass with small, but positive, specific angular momentum (yielding a maximum bulge mass) . 
As a simple intermediate case, consider the scenario in which the zero angular momentum 
bulge forms out of the material with AMD 



while the angular momentum distribution of the disk material is given by Pdisk(0 = P(0 ~ 
p{—l) where I > 0. Under the standard assumption that this material conserves its specific 
angular momentum, and that the dark matter follows an NFW density profile (Navarro, 
Frenk & White 1997), the surface density S(i?) of the resulting disk component can be 
computed using the relation 




(10) 



If the remaining material builds a disk, the resulting bulge-to-disk mass ratio is 




(11) 




(12) 



where I and R are related according to 




(13) 
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with X — R/Ryir, Vc{R) the circular velocity at radius R, and c the NFW concentration 
parameter for which we adopt c = 10 throughout. For simplicity we ignore the self-gravity 
of the disk, so that Vc{R) reflects the circular velocity of the NFW halo. 

In Figure 10 we plot the resulting disk surface density distributions S(-R) for four rep- 
resentative halos from our sample of 23 (thick solid lines). In each panel we indicate the 
bulge-to-disk ratio computed with equation (11). The thin lines correspond to disks that 
would form out of an angular momentum distribution given by 

CfiQ, 1) 



which corresponds to the 'Universal' profile of equation 9. Here 

/i 



l-{pL- l)ln 



(15) 



where Agas is the spin parameter of the gas in the halo under consideration (listed in Table 1). 
We adopt // = 1.25, which corresponds to the median value of all halos analyzed by BOl. 

Whereas the disks forming out of the AMDs of the gas analyzed here have surface 
brightness profiles that are very close to exponential (corresponding to a straight line in 
the panels on the left), this is not the case for disks forming out of the 'Universal' AMD. 
The largest difference between the two surface brightness profiles is apparent at small radii. 
This is due to the bulge-formation scenario included in one case but not the other. In 
fact, as shown by van den Bosch (2001), if a simple bulge formation scenario is included, 
the 'Universal' AMD also yields near-exponential disks. One problem that is not solved by 
turning low angular momentum material into a bulge, however, is the fact that the 'Universal' 
AMD yields disks with a clear truncation radius which occurs at too high surface density 
(van den Bosch 2001). However, the disks forming out of the AMDs of the gas analyzed here 
continue as exponentials to much lower surface brightness, suggesting that the cell averaging 
adopted by BOl underestimated the true jmax- 

Thus, the angular momentum distributions of the gas in proto-galaxies seem consistent 
with the observed surface density distributions of disk galaxies if the material with negative 
specific angular momentum is assumed to make a bulge component. Nevertheless, several 
problems remain. 

First of all, since 0.05 <, /gas ^ 0.5 in the scenario considered above, a large fraction of 
the halos will form bulge dominated 'early-type' systems. If we use a bulge-to-disk ratio of 
B/D — 1 to discriminate between late and early-type galaxies, disk dominated late type 
galaxies only form in ~ 60 percent of all halos in our sample. Furthermore, with a minimum 
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/gas of about 5 percent (see Figure 2), no systems will form with B/D < 0.1, so the question 
of how bulge-less dwarf and LSB galaxies form (van den Bosch 2001; van den Bosch et 
al. 2001) is unanswered. 

Another problem concerns the various misalignments of angular momentum. As we 
have shown in the lower panels of Figure 7, not all gas has the same direction of angular 
momentum. When the self-gravity of the gas becomes important, gas parcels with different 
angular momentum directions (moving on non-coplanar circular orbits) will exert torques 
on each other, resulting in a redistribution of angular momentum. In addition, when the 
gas is misaligned with the potential of the dark matter halo^, transfer of angular momentum 
between the gas and dark matter will occur. If the gas predominantly transfers angular 
momentum to the dark matter, which one might naively expect since the latter is the dom- 
inant mass component, this results in an additional problem as many authors in the past 
have pointed out: In order to explain the observed distribution of disk scale lengths the gas 
cannot lose much of its angular momentum. 



5. Conclusions 

It is generally assumed that the gas and the dark matter in a proto-galaxy have the 
same distribution of specific angular momentum. This is motivated by the fact that both 
mass components experience the same tidal torques. However, the two components undergo 
rather different relaxation mechanisms: whereas the dark matter undergoes collisionless, 
violent relaxation, the gas gets shocked which could in principle cause a redistribution of 
its angular momentum distribution. In this paper we have used numerical simulations to 
test, for the first time, whether the ansatz that gas and dark matter have similar initial 
angular momentum distributions is correct. We presented a numerical simulation of structure 
formation in a ACDM cosmology including both dark matter and a non-radiative gas. For 
~ 380 halos with virial masses in the range 4.8 x l^fh"'^ Mq < M^ir < 1.1 x lO^^/i"^ Mq (at 
z = 3) we computed the angular momenta of both the gas and the dark matter. 

We have shown that on average the gas and dark matter have the same distribution of 
spin parameters and that their detailed angular momentum distributions in individual halos 
are remarkably similar. Evidently, shocking during the virialization process does not decouple 



^We emphasize here that any possible misahgnment between the gas and the moment of inertia of the 
dark matter halo is not related to the misalignment angle 9 between the angular momentum vectors of the 
gas and dark matter: halos are not rotationally supported systems, and as shown by Porciani, Dekel & 
Hoffman (2002), the spin and inertia of dark matter halos are poorly aligned. 
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the angular momentum of the gas from that of the dark matter, supporting the standard 
assumption that gas and dark matter have identical angular momentum distributions. We 
have shown that these results arc robust by also analyzing a high resolution numerical 
simulation of a single cluster-sized halo at ^ = in a much larger cosmological volume. 

A careful analysis of the detailed angular momentum distributions reveals that they 
arc inconsistent with the 'Universal' profile suggested by BOl. In the angular momentum 
distributions presented here between 5 and 50 percent of the mass has negative specific 
angular momentum, whereas BOl have (implicitly) assumed that all mass has positive specific 
angular momentum. Furthermore, our angular momentum distributions have wings out to 
very high specific angular momentum, whereas the profile suggested by BOl has a cut-off 
at a certain jmax that results in disk truncation at too high surface density (van den Bosch 
2001). 

The fact that relatively large mass fractions have negative specific angular momentum 
suggests that the standard assumption of detailed specific angular momentum conservation, 
which is often made in theories of disk formation, cannot be correct. After all, disk galax- 
ies do not typically contain material with negative specific angular momentum. We have 
suggested that during the cooling process the gas with negative specific angular momentum 
may collide/shock with material with positive specific angular momentum to build a bulge 
component with zero angular momentum. In such a picture at least the total amount of an- 
gular momentum of the gas is conserved, and we have shown that it can produce disks with 
near-exponential surface density distributions without the small truncation radii implied by 
the 'Universal' profile of BOl. Numerical simulations of disk formation that include cooling 
do indeed seem to reveal the formation of a dense central knot of gas in addition to the 
disk (e.g., Katz & Gunn 1991; Navarro & White 1994). Although this 'bulge' is typically 
interpreted as forming out of sub-clumps that lose their specific angular momentum to the 
dark matter through dynamical friction, our results suggest that even without substructure 
a bulge may form naturally. In fact, taking the results at face value suggests that ~ 40 
percent of all halos are unable to form disk dominated galaxies, simply because they have 
too large mass fractions with negative specific angular momentum. In addition, systems 
with bulge-to-disk mass ratios below 10 percent will be extremely rare. 

In reality the formation of disk galaxies will be more complicated as depicted above. 
We have shown that the angular momentum vectors of the gas and dark matter are mis- 
aligned by on average ~ 35°. Surprisingly the amount of misalignment is not a measure for 
the degree to which the respective angular momentum distributions agree or disagree with 
each other. There is a weak tendency for the misahgnment to be weaker in halos that are 
more massive or have larger spin parameters. Higher resolution simulations are required 



-18- 



to investigate the amount and origin of this misahgnment in more detail. In addition to 
this misahgnment between the two components, the angular momentum vector of each in- 
dividual component can change direction quite dramatically with radius. This is typically 
interpreted as originating from dccoherence in the angular momentum orientation of mate- 
rial accreted at different times (Ryden 1988; Quinn & Binney 1992) and is often invoked as 
a mechanism to create warps in disk galaxies (e.g., Ostriker & Binney 1989). Although the 
exact implications of these various misalignments are currently unclear, and require a more 
detailed investigation, it is clear that any misalignment of the gas with either itself or with 
the inertia moment of the dark matter halo has to vanish if one wants to build a stable disk. 
If this results in a significant transfer of angular momentum from the gas to the dark matter 
these misalignments might be an additional cause for the angular momentum catastrophe 
hampering the formation of sufficiently extended disk galaxies in numerical simulations. 

Finally, we emphasize that the implications for disk/bulge formation discussed above 
are highly speculative. In our simulations the gas was not allowed to cool, and we used 
the angular momentum distributions of the gas found at z = 3 to speculate about the fate 
of the gas once cooling were turned on. In reality the gas will already have (partially) 
cooled inside sub-clumps that merge to form the systems analyzed here, and this will also 
impact on the angular momentum distribution of the gas and its subsequent evolution (e.g., 
dynamical friction). On the other hand, it has been argued by several authors that in order 
to prevent the angular momentum catastrophe the gas has to be prevented from cooling 
inside sub-clumps for instance by feedback processes or preheating (e.g., Navarro & White 
1994; Weil, Eke & Efstathiou 1998; Sommer-Larsen, Gelato & Vedel 1999). Under those 
circumstances the angular momentum distribution of the gas may in fact resemble those 
analyzed here, although it is unclear to what extent these heating processes influence the 
angular momentum distributions of the gas prior to being incorporated in the halos. For 
instance, processes related to feedback and re-ionization create pressure gradients in the gas 
which may decouple the angular momentum distribution of the gas from that of the dark 
matter. In a follow-up paper (Abel et al. 2002) we address in more detail the impact such 
pressure forces may have on the angular momentum distributions in proto-galaxies. 
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to James Bullock, Andreas Burkert, Avishai Dekel, Ariyeh Mailer, Joel Primack and Risa 
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Fig. 1. — The distributions of Adm (right panel) and Agas (left panel) for all 378 halos in our 
sample. Solid lines are the best fit log-normal distributions of equation (6). The values of 
the best fit parameters are indicated in each panel. Note that the gas and dark matter have 
very similar distributions of spin parameters. 
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Fig. 2. — A comparison of various global angular momentum parameters of the dark matter 
and gas for all 378 halos in our sample. The upper left panel plots Agas versus Adm- The 
dotted line corresponds to perfect agreement (i.e., Agas = Adm) and is plotted for comparison. 
Although on average the spin parameters of gas and dark matter are similar, for each indi- 
vidual halo the difference between the spin parameters of both components can be relatively 
large. In the upper right panel we plot f^^^ versus /dm- Again, on average the mass fractions 
with negative specific angular momentum are similar for the gas and the dark matter (i.e., 
(/oM/Zgas) — 0.97 ± 0.10). The two panels in the middle plot f^^ versus Agas (left) and /^j^ 
versus Adm (right). Note how the total specific angular momentum of each component is 
anti-correlated with its respective mass fraction with negative jz- The same is true when 
plotting /gas rather than f^^ (lower left panel). Finally, the lower right panel plots /gas versus 
/gag. As expected, all halos have /gas < /gas- Note however that the difference depends on 
/gag and thus on Agas- See the text for a more detailed discussion. 



-24- 




Fig. 3. — The distribution of the angle 9 between the total angular momentum vectors of 
the dark matter and the gas. The average and mean of the distribution are indicated. 
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Fig. 4. — The upper left panel plots Adm — ^gus as function of the halo virial mass. The thick 
solid line in the middle indicates the run of the average Adm — Agas, while the upper and lower 
solid lines indicate the 68 percent interval of Adm — Agas at fixed M^i^.. Note that the difference 
in the spin parameters of the gas and the dark matter do not depend significantly on the 
halo virial mass. The other panels plot 9 as function of the absolute difference between the 
spin parameters of gas and dark matter (upper right panel) , the total virial mass (lower left 
panel), and the spin parameter of the gas (lower right panel). Thick solid lines correspond to 
the mean and 68 percent interval of 9. Whereas there is no dependence of 9 on |Adm — Agas|, 
the average misalignment between the angular momentum vectors of the gas and dark matter 
is larger in less massive halos and in halos with smaller spin parameters. 
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l''=jV(RvirVvir) 

Fig. 5. — The angular momentum distributions p{l'") for the 20 halos hsted in Table 1. Sohd 
and dashed lines correspond to the dark matter and gas, respectively. Note the remarkable 
similarity between the angular momentum distributions of the gas and the dark matter. 
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l=j/(RvlrVvlr) 

Fig. 6. — A comparison of the angular momentum distributions p{l) (solid lines) and p{l^) 
(dashed lines) of the gas in the 20 halos listed in Table 1. Note how p{l) typically has much 
smaller wings at negative specific angular momentum than p(/^) (cf, lower right panel of 
Figure 2). 
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Fig. 7. — The values of the spin parameter inside radius R as functions of R/Rv\r for the 20 
halos hsted in Table 1. Solid (dashed) lines correspond to the dark matter (gas). Typically 
\{R) either decreases with radius or stays roughly constant. Although in some halos Adm(-R) 
and Agas(i?) can be quite different, there is no indication for a systematic difference between 
\{R) of the gas and dark matter. 
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Fig. 8. — The misalignment angle A9{R) between the angular momentum vector inside 
radius R and that of the entire halo. As in Figure 7 sohd (dashed) hues correspond to 
the dark matter (gas). Results are plotted as function of i?/i?vir- Note that the angular 
momentum directions of both the gas and dark matter can vary quite significantly with 
radius, and in some cases is clearly different for the gas and the dark matter. 
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Fig. 9. — Results at 2; = for the high resolution simulation of cluster 'SO' described in 
Section 3.3. The upper left panel plots for both the dark matter (solid line) and the 
gas (dashed line). As for the halos presented in Figure 5 the dark matter and gas have 
virtually indistinguishable distributions of specific angular momentum. The upper right 
panel compares p(/) (solid line) with piT) for the gas (cf. Figure 6). Note that even when 
considering the streaming motions of the gas, there is a large fraction of gas mass with 
negative specific angular momentum (see also Table 1). The lower two panels plot A(i?) and 
A^(i?) for halo 'SO', where the sofid (dashed) lines correspond to the dark matter (gas). A 
comparison with Figures 7 and 8 shows that halo 'SO' behaves similar as the halos analyzed 
in Section 3.2. 
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Fig. 10. — The thick sohd hnes correspond to the density distributions (in arbitrary units) 
of the disks that form in four halos hsted in Table 1 under the assumptions that (i) material 
with negative specific angular momentum forms a zero-angular momentum bulge together 
with material with the same, but positive, specific angular momentum distribution (e.g., 
equation [10]), (ii) the remaining material forms a disk under detailed angular momentum 
conservation, (iii) the dark matter halo has an NFW density distribution with c = 10, and 
(iv) the self-gravity of the disk can be ignored. The resulting bulge-to-disk mass ratios 
B/D are indicated in each panel. The thin solid lines correspond to disks that form out 
of an angular momentum distribution given by the 'Universal' profile of equation (14) with 
/i = 1.25 and Agas as listed in Table 1. We plot Sdisk(-R) both as function of R/R^ij. (panels 
on the left) and \og{R/ R^ij.) (panels on the right), in order to emphasize the differences at 
large and small radii, respectively. 
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Table 1: Glo bal Parameters of Halos in Figures 5-8. 
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Note. — Global parameters of the halos whose angular momentum distributions are presented in Figures 5- 
8. Column (1) lists the halo ID. Columns (2) and (3) list the number of dark matter and gas particles, 
Columns (4) and (5) the spin parameters, and Columns (6) and (7) the mass fractions with negative j^. 
For comparison, column (8) lists the gas mass fraction with negative jz ■ Finally, column (9) lists the 
misalignment angle 9. 



